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Over the past few years it has been discovered that an "observable" can be 
set up on the lattice which obeys the discrete Cauchy-Ricmann equations. 
The ensuing condition of discrete holomorphicity leads to a system of linear 
equations which can be solved to yield the Boltzmann weights of the underlying 
lattice model. Surprisingly, these are the well known Boltzmann weights which 
satisfy the star-triangle or Yang-Baxter equations at criticality. This connection 
has been observed for a number of exactly solved models. I briefly review these 
developments and discuss how this connection can be made explicit in the 
context of the Zjv model. I also discuss how discrete holomorphicity has been 
used in recent breakthroughs in the rigorous proof of some key results in the 
theory of planar self-avoiding walks. 
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1. Introduction 

The study of exactly solved lattice models in statistical mechanics has been 
a richly rewarding area of activity over many decades. Among the pan- 
theon of exactly solved models is the two-dimensional Ising model, which is 
fundamental to the theory of phase transitions and critical phenomena. 1-3 
Two-dimensional lattice models with local interactions like the Ising model 
possess remarkable properties at criticality. The fluctuations in such sys- 
tems are conformally invariant at the critical point. One of the triumphs 
of statistical physics is that this property entails the classification of uni- 
versality classes of second order phase transitions in terms of conformal 
field theories. 3 Roughly speaking, there is an exactly solved model for each 
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universality class. The beauty of exactly solved models is that their critical 
Boltzmann weights satisfy the star-triangle or Yang-Baxter equation. 1-3 

A remarkable link has been uncovered between Yang-Baxter integra- 
bility and the condition of discrete holomorphicity. 4-6 This involves the 
establishment of a local identity which is satisfied by an operator defined 
on the lattice. This operator is an example of a discrete parafcrmionic ob- 
servable appearing in conformal field theory. The observable obeys a lattice 
version of the Cauchy-Riemann equations and is thus known to be dis- 
cretely holomorphic. The condition of discrete holomorphicity leads to a 
system of linear equations which can be solved to yield the Boltzmann 
weights of the lattice model. These are the Boltzmann weights which are 
known to satisfy the Yang-Baxter equations at criticality. This connection 
has been observed for a number of exactly solved models, 4-6 including the 
Zat model, 4 the Potts model 5 and the O(n) loop model. 5 

Here we will discuss these developments in terms of the O(n) loop model. 



2. The O(n) loop model 

The O(n) loop model on the square lattice was originally constructed by 
Nienhuis. 7 ' 8 It is exactly solved by means of the Bethe Ansatz, with pe- 
riodic 9 and open 10,11 boundaries. As an integrable model, its underlying 
i?-matrix is that of the Izergin-Korepin, 12 or vertex model. 10 The par- 
tition function is given by 

9 

which is a sum over all possible loop configurations Q. Here rrij is the number 
of occurrences of Boltzmann weight pi and P is the total number of closed 
loops of fugacity n in any given configuration. The nine non-zero Boltzmann 
weights pi are shown in Fig. 1. These weights are given by 7,10 

pi = sin(3A — u) sin u + sin 2A sin 3A 
p2 = P3 = £i sin(3A — u) sin2A 
Pi = P5 = C2sinusin2A 

P6 = Pi = sin(3A — u) sin it (2) 
pa = sin(3A — u) sin(2A — u) 
p g = — sin(A — u) sin u 

with e\ = e% = 1, < u < 3A and n = -2cos4A. 
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Fig. 1. The Boltzmann weights of the O(n) loop model. 

The parafermion operator is defined on the lattice as 5 

F(z) = P(G)e-' laW{z) . (3) 

eer(o,z) 

Here P{G) is the probability of configuration Q, r(0, z) is a set of loop 
configurations for which points and z belong to the same loop, W(z) is 
the winding angle of this loop from points and z, and a is the spin of the 
parafermion F(z). For the 0(n) loop model 

F(z)= e-^^^Vr ■■■pTn P (4) 

7(a— «)6f! 

for some path 7(0 — > z) across a domain J7. 

What Ikhlef and Cardy 5 have shown is that the parafermion F(z) is 
discretely holomorphic, i.e., it satisfies 

F(p) - e w F{q) - F(r) + e w F(s) = (5) 

on each vertex as defined in Fig. 2. This condition is equivalent to requiring 
that F is divergence and curl free. 




Fig. 2. The labelling of a vertex. 

To establish this, one chooses a particular mid edge p. There are four 
possible different external connectivities, as shown in Fig. 3. These possi- 
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bilities result in four equations 

Pi + Vp-2 - V(,~ X P4 - P7 = 

-C _1 P2 + np 5 + v(p 7 - ^C 1 (ps + n p 9 ) = 

n P3 - CPi - v(~ 2 P7 + v(nps + pg) = (6) 
-K~ 2 P2 + ^CP4 + np 6 - C 2 ps - C 2 P9 = 

where v = e l6 '^' T+1 ' and C = e ICT7r . The remarkable point is that Ikhlef and 
Cardy 5 were able to solve this set of linear equations to yield the Boltzmann 
weights px,...,pg and the conformal spin 



provided u = a (6 — n) + 9. In particular, at 9 = ^, u = 3A/2, the isotropic 
point. Substitution of the Boltzmann weights (2) into equations (6) shows 
that these equations are indeed satisfied. 




Fig. 3. The four possible external connectivities from mid edge p. 

The surprising point is that the Yang-Baxter equations are cubic in the 
Boltzmann weights, whereas imposing discrete holomorphicity led to a set 
of equations which are linear in the Boltzmann weights. There thus appears 
to be a simpler route to uncovering integrable lattice models at criticality. 
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Why is this? And what is the precise link between discrete holomorphicity 
and integrability? These questions will be taken up further in Section 4 
in the context of the Zjv model, which can be discussed from an algebraic 
point of view. 

We now turn to the related developments in the theory of self-avoiding 
random walks. 

3. Self-avoiding random walks 

The random walk, also known as the drunkard's walk, was introduced over 
100 years ago. 13 The study of random walks has since developed into the 
central pillar of the theory of probability and random processes. 14 A key 
feature of the random walk is that, like the drunkard, it has no memory of 
where it has been. This absence of memory is the defining property of a 
Markovian process. 

A more challenging class of mathematical problems arise when a random 
walk is not allowed to revisit a site it has already visited. Such self-avoiding 
random walks were originally proposed as models of long polymer chains in 
a good solvent. Despite the simplicity of their definition, the non-Markovian 
nature of self-avoiding random walks makes them notoriously difficult to 
study. Nevertheless much is known about the theory of self-avoiding walks, 
particularly in higher dimensions. 15 Much is also known in two and three 
dimensions from extensive computer enumeration. The self-avoiding con- 
straint is more restrictive in two dimensions and thus self-avoiding walks 
on planar lattices are arguably more interesting to study. 




(a) (b) (c) 

Fig. 4. The (a) honeycomb, (b) 3-12 and (c) martini lattices. 

In 1982 Nienhuis 16 was able to obtain some key information for self- 
avoiding walks by considering the O(n) loop model on the honeycomb 
(hexagonal) lattice depicted in Fig. 4(a). In general, the O(n) model has a 
connection to self-avoiding walks in the limit n — 0. 17 The critical Boltz- 
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mann weights of the model on the honeycomb lattice were obtained from an 
exact mapping to the Potts model. 16 The Boltzmann weights defining the 
critical line of the O(n) loop model on the honeycomb lattice also follow 
as a condition for the model to be solved exactly in terms of the Bethe 
Ansatz. 18-20 More generally the O(n) loop model on the honeycomb lattice 
follows as a limiting case of the exactly solved O(n) loop model on the 
square latticed 

The critical Boltzmann weights of the model give the value for the con- 
nective constant /i, which is a measure of the growth in the number of 
allowed self- avoiding walks with increasing length. The Nienhuis result, 



2 + s/2 = 1.847 759 . . . (8) 



was generally believed to be exact, having been confirmed to high precision 
from extensive enumeration of self-avoiding walks on the honeycomb lattice. 
No such exact results are known for the other regular planar lattices. 

Similar exact results are known for self-avoiding walks on related semi- 
regular lattices. These are the 3-12 lattice (Fig. 4(b)) and the martini lattice 
(Fig. 4(c)). For the 3-12 lattice, the exact value 21 ' 22 pt = 1.711 041 .. . follows 
from the solution of 

111 (9) 



M M V2 + V2 

And for the martini lattice, the exact value 23 /i = 1.750 564 . . . follows from 
the solution of 

\ + -,= , 1 n - (io) 

M M V2 + V2 

More generally, the 0(n) model on the honeycomb lattice can be mapped 
to the O(n) model on the 3-12 22 and martini lattices. 23 

Up until recently, a rigorous mathematical proof of the result (8) had 
seemed out of reach. In a remarkable paper, Duminil-Copin and Smirnov 24 
have given an ingenious proof of the Nienhuis result. The first part of the 
proof involves defining a parafermionic observable on the honeycomb lattice 
which is shown to be discretely holomorphic. The Boltzmann weights of the 
model follow from the condition of discretely holomorphicity. In this sense 
the value (8) also appears as a consequence of discrete holomorphicity. The 
second part of the proof builds on previously hard won results in the general 



a Specifically, when either u = A or u = 2A in (2), corresponding to the two orientations 
in Fig. 5. 
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theory of self-avoiding walks to rigorously establish that this value is indeed 
the connective constant. 



3.1. The polymer adsorption transition 

These arguments have been extended 25,26 to provide a rigorous proof of the 
critical surface fugacity y c of self-avoiding walks on the honeycomb lattice 
in the presence of a boundary. This is the critical adsorption temperature 
at which the polymer becomes attached to the surface. For y > y c the 
polymer is adsorbed and for y < y c the polymer is desorbed. The critical 
value depends on the orientation of the lattice (see Fig. 5). 




(a) (b) 



Fig. 5. The two orientations of the honeycomb lattice. 



For the orientation shown in Fig. 2(a), the exact value 27 

y c = l + V2 = 2.414... (11) 

had been obtained earlier using the underlying integrability of the O(n) 
loop model with a boundary. 10 For the orientation shown in Fig. 2(b), the 
exact value 28 



y c = J 2 + >/2 = 

V 1 + y/2 - V 2 + V2 
also followed from the integrable boundary weights of the O(n) loop 
model. 10 

Building on the approach of Duminil-Copin and Smirnov, these results 
have been shown by Beaton et al. 25,26 to also follow from the discrete holo- 
morphicity of a parafermionic observable. Their mathematical proof that 
these values are indeed the critical adsorption temperatures is an impres- 
sive tour de force, involving the adaption of previous rigorous work on 
self-avoiding walks and the rigorous treatment of self-avoiding bridges in a 
finite strip. 

We turn now to further examine the connection between discrete holo- 
morphicity and integrability. 
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4. The Zjv model 

The Z^v model is an TV-state model in which the spins s r at vertex r take 
values from the Nth roots of unity, s r = uj qr , with q r e {0,1,..., N— 1} and 
lo = exp(27ri/A^). The model contains a number of well known models as 
special cases, namely the Ising model (N = 2), the three-state Potts model 
(N = 3) and the Ashkin- Teller model (N = 4). The Z N model has a well 
known set of critical points obtained by Fateev and Zamolodchikov (FZ) as 
solutions to the star-triangle relations. 29 By defining lattice parafermions 
as products of neighbouring order and disorder variables with a suitable 
phase factor, it was shown by Rajapour and Cardy that the FZ Boltzmann 
weights ensure that the lattice parafermions obey the discrete version of 
the Cauchy-Riemann equations, i.e., they enjoy the property of discrete 
holomorphicity. 4 ' 6 These results were extended to lattices with anisotropic 
interactions. 

The underlying motivation for this work on the Zjy model is the con- 
struction of discretely holomorphic parafermions, which are expected to 
become the holomorphic parafermions of the conformal theory. Here our in- 
terest is in the connection with integrability. From that perspective, the key 
finding is that the contour sum around each elementary face or rhombi of the 
lattice vanishes, the key ingredient for which is the set of FZ critical Boltz- 
mann weights. This approach was extended by explicitly considering the 
condition of discrete holomorphicity on two and three adjacent rhombi. 30 
For two rhombi this leads to a quadratic equation for the Boltzmann weights 
which implies the known inversion relations for the Zjy model. For three 
rhombi it leads to a cubic equation from which the star-triangle relation 
follows. The simplicity of the discrete holomorphic approach is that the 
two- and the three-rhombus equations are equivalent to the one-rhombus 
equation. 

5. Discussion 

As alluded to in the previous section, the original motivation behind the 
recent progress on various aspects of discrete holomorphicity covered in 
this article was to set up the necessary discrete building blocks for the 
passage to the continuous scaling limit. There is ample motivation for 
doing this. For example, the geometry of the curves traced out by self- 
avoiding walks is conjectured to be described by the mathematical theory 
of Schram-Loewner evolution (SLE 8 / 3 ). 31 Basically, if the scaling limit of 
two-dimensional self-avoiding walks exists and has a certain conformal in- 
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variance property, then the scaling limit must be described by the particular 
stochastic process SLE 8 / 3 . However, establishing the existence of the scal- 
ing limit (and therefore also its conformal invariance) is a difficult open 
problem. Nevertheless, significant progress in this direction has been made 
for other two-dimensional models, most notably percolation and the Ising 
model. 32 - 34 

Although the programme to establish the scaling limit of self-avoiding 
walks in a mathematically rigorous way is not yet complete, the first step, 
involving establishing analytic properties on the lattice via the discrete 
Cauchy-Riemann equations, has led to some remarkable developments. 
These include the surprising mathematical proofs of the results (8), (11) 
and (12) for the connective constant and critical surface fugacity for self- 
avoiding walks on the honeycomb lattice. 24-26 

There are also ongoing developments with regard to the connection with 
integrability, particularly at a boundary. A boundary condition can also be 
imposed on a discretely parafermionic observable. 35 ' 36 This approach has 
led to a set of new integrable boundary weights for the Zjy model. 35 Several 
new sets of boundary weights have also been obtained in this way for the 
O(n) and loop models. 36 

More developments are eagerly awaited. 
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